V odd fluctuations and Long Range Order in Heavy Ion Collisions. 
Deformed QCD as a Toy Model. 
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We study the local violation of T 3 and CP invariance in heavy ion collisions as observed at RHIC 
and LHC using a simple "deformed QCD" model. This model is a weakly coupled gauge theory, 
which however has all the relevant crucial elements allowing us to study difficult and nontrivial 
questions which are known to be present in real strongly coupled QCD. Essentially, we want to un- 
derstand the physics of long range order in form of coherent low dimensional vacuum configurations 
observed in Monte Carlo lattice simulations. Apparently precisely such kind of configurations are 
responsible for sufficiently strong intensity of asymmetries observed in heavy ion collisions. 



I. INTRODUCTION AND MOTIVATION 

Recently it has become clear that quantum anomalies 
play very important role in the macroscopic dynamics of 
relativistic fluids. Much of this progress is motivated by 
very interesting ongoing experiments on local V and CP 
violation in QCD as studied at RHIC and ALICE at the 
LHC [1-5]. It is likely that the observed asymmetry is 
due to charge separation effect [6, 7] as a result of the 
chiral anomaly, though some background processes may 
also contribute to the measured observed asymmetry [8- 
12]. The ideas formulated in [6, 7] were further developed 
in follow up papers [13, 14] where the effect was coined 
as chiral magnetic effect (CME). 

We shall not discuss a number of subtle questions of the 
CME in the present work by referring to a recent review 
[15]. Instead, we concentrate on a single crucial element 
for CME to be operational. Namely, the key assump- 
tion of the proposal [6, 7] is that the region where the 
so-called {9(x,t) in d} ^ should be much larger in size 
than the scale of conventional QCD fluctuations which 
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have typical correlation lengths of order -~ ^^q CD - 
9(x, t)i n d parameter enters the effective lagrangian as fol- 
lows, C e = -6 md q where q = ^ze flI/pa G afiu G a( " T is the 
topological density operator, such that local V and CV 
invariance of QCD is broken on the scales where corre- 
lated state with (9(x, t)i n d) ^ is induced. As a result of 
this violation, one should expect a number of V and CV 
violating effects taking place in a relatively large region 
where (9(x,t)i n d) ^ 0. The question which is addressed 
in the present work is as follows: what is the physics 
behind of this long range order when {6(x, t)i n d) is corre- 
lated on distances much larger than conventional Aq CD 
scale? 

One should say that such kind of long range order 
indeed has been observed in recent lattice simulations 
which have revealed some very unusual features. To be 
more specific, the gauge configurations studied in [16-19] 
display a laminar structure in the vacuum consisting of 
extended, thin, coherent, locally low-dimensional sheets 
of topological charge embedded in 4d space, with oppo- 
site sign sheets interleaved. A similar structure has been 
also observed in QCD by different groups [20-25] and also 



in two dimensional CP N 1 model [26]. 

The key element relevant for the present studies is 
the observation that the important gauge configurations 
which contribute to the CME are in fact very similar 
to extended, locally low-dimensional sheets of topologi- 
cal charge embedded in 4d space studied in refs. [16-25]. 
Such a similarity was noticed previously using numerical 
lattice simulations [24, 25] . The main goal of the present 
work is to understand this relation on a deeper theoret- 
ical level using "deformed QCD" as a toy model. Fur- 
thermore, we shall argue that an external magnetic field 
(which is a key element in studying CME) may serve as 
a complimentary tool to analyze the relevant topological 
configurations present in the QCD vacuum. 

We start in section II by reviewing the relevant parts 
of the model [27, 28]. This is a simplified version of QCD 
which, on one hand, is a weakly coupled gauge theory 
wherein computations can be performed in theoretically 
controllable manner. On other hand, the correspond- 
ing deformation preserves all the relevant elements of 
strongly coupled QCD such as confinement, degeneracy 
of topological sectors, nontrivial 9 dependence, and many 
other important aspects which allow us to test some fas- 
cinating features of strongly interacting QCD, including 
the CME and long range order as we shall argue below. 

In section III we introduce the matter field into the 
model such that it couples to the external £7(1) magnetic 
field. We observe that the matter field essentially traces 
the long range order which is inevitable feature of pure 
gauge theory. 

In section IV we apply these ideas to CME. We observe 
that CME takes place precisely on the same extended 
configurations which arc identified with low-dimensional 
sheets studied in lattice simulations [16-25]. In different 
words, the CME which takes place in the presence of the 
magnetic field precisely traces the structure of the orig- 
inal topological configurations which are always present 
in QCD vacuum. This observation suggests that an ex- 
ternal magnetic field could be a perfect complimentary 
tool which allows to study a very complicated QCD vac- 
uum structure which apparently shows a long range order 
rather than some localized finite size fluctuations such as 
instantons. 

Section V is our conclusion where we consider appli- 
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cations of these ideas to V odd fluctuations observed at 
RHIC and ALICE at the LHC [1-5]. We interpret a rela- 
tively strong asymmetry observed in heavy ion collisions 
as a direct manifestation of the long range order studied 
in present work. We argue that the asymmetry would be 
much weaker than observed if this long range order was 
not present in the system. 



Parameter L here is the length of the compactified dimen- 
sion which is assumed to be small. The coefficients of the 
polynomial P [f2(x)] can be suitably chosen such that the 
deformation potential (2) forces unbroken symmetry at 
any compactification scales. At small compactification L 
the gauge coupling is small so that the scmiclassical com- 
putations arc under complete theoretical control [27]. 



II. DEFORMED QCD 



B. Infrared description 



Here we overview the "center-stablized" deformed 
Yang-Mills developed in [27] and reviewed in II A. In the 
deformed theory an extra term is put into the Lagrangian 
in order to prevent the center symmetry breaking that 
characterizes the QCD phase transition between "con- 
fined" hadronic matter and "deconfined" quark-gluon 
plasma. Thus we have a theory which remains con- 
fined at high temperature in a weak coupling regime, 
and for which it is claimed [27] that there does not exist 
an order parameter to differentiate the low temperature 
(non-abelian) confined regime from the high temperature 
(abelian) confined regime. One should remark here that 
"deformed QCD" model has been previously successfully 
used to test some highly nontrivial features of strongly 
coupled QCD such as emergence of non-dispersive con- 
tact term in topological susceptibility [28] and emergence 
of the topological Casimir behaviour in gauge theory with 
a gap [29]. Finally, what is most relevant for present 
studies is the observation that "deformed QCD" model 
exhibits a long range order and double layer structure 
which was interpreted in ref. [30] as a manifestation of 
the corresponding features observed in strongly coupled 
regime in the lattice simulations. The extended structure 
in this model is manifested in a form of the domain walls, 
which are present in the system as a result of a generic 2-7T 
periodicity of the effective low energy Lagrangian govern- 
ing the dynamics of "deformed QCD". The low energy 
description will be reviewed in section II B while the part 
on domain walls in pure glue theory will be reviewed in 
section II C. 



A. The model 

We start with pure Yang-Mills (gluodynamics) with 
gauge group SU{N) on the manifold R 3 x S* 1 with the 
standard action 

S YM = I d* x -1 tr [**,(*)] , (1) 
and add to it a deformation action, 



AS= [ d 3 x -^P [fi(x)] 



(2) 



built out of the Wilson loop (Polyakov loop) wrapping 
the compact dimension 



fi(x) =V \e i $ d 



C4 A4(x,X4) 



(3) 



As described in [27], the proper infrared description 
of the theory is a dilute gas of N types of monopoles, 
characterized by their magnetic charges, which are pro- 
portional to the simple roots and affine root a a G A a ff 
of the Lie algebra for the gauge group U(l) N . For a fun- 
damental monopole with magnetic charge a a e A a g, the 
topological charge is given by 
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and the Yang-Mills action is given by 
d^x-} 
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The ^-parameter in the Yang-Mills action can be included 
in conventional way, 
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with F^ = e^P a F p<y . 

The system of interacting monopoles, including 9 pa- 
rameter, can be represented in the dual sine-Gordon form 
as follows [27, 28], 
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where £ is magnetic monopole fugacity which can be 
explicitly computed in this model using conventional 
scmiclassical approximation. The 9 parameter enters 
the effective Lagrangian (7) as 9/N which is the direct 
consequence of the fractional topological charges of the 
monopoles (4). Nevertheless, the theory is still 2ir peri- 
odic. This 2ir periodicity of the theory is restored not 
due to the 2ir periodicity of Lagrangian (7). Rather, it is 
restored as a result of summation over all branches of the 
theory when the levels cross at 9 = ir(mod 27r) and one 
branch replaces another and becomes the lowest energy 
state as discussed in [28]. 

Finally, the dimensional parameter which governs the 
dynamics of the problem is the Debye correlation length 
of the monopole's gas, 



ml = LC 



(8) 
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The average number of monopoles in a "Debye volume" 
is given by 



form [30] 
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(9) 



The last inequality holds since the monopole fugacity is 
exponentially suppressed, £ ~ e^ 1 ^ 9 , and in fact we can 
view (9) as a constraint on the validity of the approxi- 
mation where semiclassical approximation is justified. 



C. Domain Walls in deformed QCD 

A discrete set of degenerate vacuum states as a re- 
sult of the 2tt periodicity of the effective Lagrangian 
(7) for er field is a signal that the domain wall config- 
urations interpolating between these states are present 
in the system. However, the corresponding configura- 
tions are not conventional static domain walls similar to 
the well known ferromagnetic domain walls in condensed 
matter physics which interpolate between physically dis- 
tinct vacuum states. In contrast, in present case a cor- 
responding configuration interpolates between topologi- 
cally different but physically equivalent winding states 
| re), which are connected to each other by large gauge 
transformation operator. Just because of that, the corre- 
sponding domain wall configurations in Euclidean space 
should be interpreted as configurations describing the 
tunnelling processes in Minkowski space, similar to Eu- 
clidean monopoles which also interpolate between topo- 
logically different, but physically identical states, see de- 
tail discussions in [30] and references therein. 

One should remark that a formal similar construction 
has been considered previously in relation with the so- 
called N = 1 axion model [31-33], and more recently in 
QCD context in [34] and in high density QCD in [35]. In 
previously considered cases [34, 35] as well as in recently 
considered case [30] there is a single physical unique vac- 
uum state, and interpolation corresponds to the transi- 
tion from one to the same physical state. 

There arc TV different domain wall (DW) types. How- 
ever, there are only (N — 1) physical propagating scalars 
er in the system as one singlet scalar field, though it re- 
mains massless, completely decouples from the system, 
and does not interact with other components at all [27]. 

In what follows, without loosing any generalities, we 
consider N = 2 case. In this case there is only one phys- 
ical field x which corresponds to a single diagonal com- 
ponent from the original SU(2) gauge group. The action 
(7) becomes, 
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x(z) = 4 arctan [exp(m CT 2;) 



(11) 



where we take 8 = to simplify things. The width of 
the domain wall is obviously determined by m" 1 , while 
the domain wall tension a for profile (11) can be easily 
computed and it is given by 



a = 2 • 
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(12) 



The topological charge density for profile (11) assumes a 
double layer structure as discussed in details in [30] 



, , C , N 4C (e m ° z - e - m " z ) 
*w I L (e m ° z + e - m " z ) 2 



(13) 



From eq. (13) one can explicitly see that the net topolog- 
ical charge Q ~ J_ ao dzq(z) on the domain wall obviously 
vanishes. However, the charge density is distributed not 
uniformly. Rather, it is organized in a double layer struc- 
ture. It was conjectured in [30] that this structure is in 
fact a trace of a similar configurations studied in the lat- 
tice simulations. 

One should also comment that, formally, a similar soli- 
ton -like solution which follows from action (10) appears 
in computation of the string tension in 3d Polyakov's 
model [27, 36]. The solution considered there emerges as 
a result of insertion of the external sources in a course 
of computation of the vacuum expectation of the Wilson 
loop. In contrast, in our case, the solution (11) is in- 
ternal part of the system without any external sources. 
Furthermore, the physical meaning of these solutions are 
fundamentally different: in our case the interpretation of 
solution (11) is similar to instanton describing the tun- 
nelling processes in Minkowski space, while in computa- 
tions [27, 36] it was an auxiliary object which appears in 
the course of computation of the string tension. 

As we already mentioned the DW described here are 
classically stable objects, but quantum mechanically they 
decay as a result of tunnelling processes. The correspond- 
ing rate was estimated in [30] and is given by 



cxp 



3 ln 2 (-M 



V4tt/ ^JUl 
cxp (-7 -M In 2 M) <C 1, 



(14) 



while the equation of motion and its solution take the 



with 7 being some numerical coefficient and M 3> 1 is 
large parameter of the model (9). The estimate (14) jus- 
tifies our treatment of the domain walls as the stable 
objects when the weak coupling regime is enforced by 
the deformation (2). 

We conclude this short overview by emphasizing that 
the most important lesson from analysis [30] is that the 
double layer structure represented by cq. (13) naturally 
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emerges in construction of the domain walls in weak cou- 
pling regime in deformed QCD. As claimed in [27] the 
transition from high temperature weak coupling regime 
to low temperature strong coupling regime should be 
smooth without any phase transitions on the way. As a 
consequence of this smoothness, the double layer struc- 
ture (13) was interpreted in [30] as a trace of the dou- 
ble layer structure observed in the lattice simulations in 
strong coupling regime [16-19]. In different words, as the 
transition from weak coupling regime to strong coupling 
regime should be smooth in this model it is naturally to 
assume that the domain walls (11) become very clumpy 
with large number of folders. Such fluctuations obvi- 
ously increase the entropy of the DW which eventually 
may overcome the intrinsic tension (12). If this happens, 
the DWs would look like as very crumpled and wrinkled 
objects with large number of foldings and effectively van- 
ishing tension as a result of large entropy. In this case an 
arbitrary number of such objects can be formed and they 
can percolate through the vacuum, forming a kind of a 
vacuum condensate. Furthermore, the DWs may loose 
their natural dimensionality, and likely to be character- 
ized by a Hausdorff dimension as recent lattice simula- 
tions suggest [24]. Nevertheless, the topological charge 
distribution with such striking features as extended co- 
herence along x, y directions and double layer structuring 
along z direction, see eq.(13), should still persist as the 
transition from weak to strong coupling regime should be 
smooth. It is quite likely that an appropriate description 
for this physics should be formulated in terms of holo- 
graphic dual model as argued in [37], however we leave 
this subject for future studies. 



combination 9^9 — rf '. Therefore we have, 



III. DOMAIN WALLS IN THE PRESENCE OF 
MATTER FIELD 



The ultimate goal of the present work is to understand 
the long range structure described above in section II C in 
the presence of external magnetic field, which is precisely 
the environment relevant for study of the CME. However, 
the gluons in pure glue theory represented in this model 
by low energy effective Lagrangian (7) do not couple to 
physical magnetic field. Therefore, we introduce a single 
charged masslcss quark field t/j into the model. The low 
energy description of the system in confined phase with a 
single quark is accomplished by introducing the rf colour 
singlet, charge neutral, meson. As usual, the 77' would 
be conventional massless Goldstone boson if the chiral 
anomaly is ignored. In the dual sine- Gordon theory the 
77' field appears exclusively in combination with the 9 
parameter as 9 — > 9 — rf as a consequence of the Ward 
Identities. Indeed, the transformation properties of the 
path integral measure under the chiral transformations 
tp — > exp(i75^-)7/> dictate that 77' appears only in the 
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where dimensionless numerical coefficient c ~ 1 can be, 
in principle, computed in this model, though our results 
do not depend on its numerical value. One can explicitly 
compute the topological susceptibility \ m this model, 
and check that the Ward Identities are automatically sat- 
isfied when the 77' field enters the Lagrangian precisely in 
form (15), see the detail computations in [28]. The 77' 
mass computed from (15) has an extra 1/N suppression, 
as it should, in comparison with mass of the er field 
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As we already mentioned, the system (15) is 2ir pe- 
riodic with respect 9 — > 9 + 2ir. This 2ir periodicity of 
the theory is not explicit in eq. (15), but nevertheless 
is restored as a result of summation over all branches of 
the theory as explained in [28]. The fact that 77' field 
enters the low energy effective Lagrangian precisely in 
combination (9 — 77') implies that system (15) is also pe- 
riodic with respect to shift rf — > rf + 2tt when summa- 
tion over all branches of the theory is properly imple- 
mented. In different words, the system (15) supports 
77' (2) domain walls as a result of a very generic feature 
of the theory 1 . This argument has been previously used 
to construct the 77' domain walls in context of the ax- 
ion physics [34]. One can derive the equation of mo- 
tion for set of rf[z) and cr(z) fields determined by La- 
grangian (15) and impose an appropriate boundary con- 
ditions rf (z = +00) — rf ' (z = -co) = 2ix to analyze this 
system numerically, in close analogy with procedure used, 
though in a different context, in ref. [34]. 

However, for our present work it is sufficient to quali- 
tatively describe the behaviour of the system in the limit 
when rn^, /m 2 G <C 1, when the basic features can be easily 
understood even without numerical computations. The 
most important lesson of this qualitative analysis, as we 
shall see in a moment, is that the light 77' field traces the 
<t(z) field by exhibiting a similar double layer structure 
discussed for pure gauge theory in section II C. 

The low energy Lagrangian which describes the lightest 
degrees of freedom for SU (2) gauge group is governed by 



1 There are many different types of the domain walls which are 
supported by the Lagrangian (15). We leave this problem of 
classification of the DWs for a future study. In this work we con- 
centrate on a simplest possible case with N = 2 to demonstrate 
few generic features of the system. 
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the following action 
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where we inserted the rj' field into eq. (10) exactly in the 
form consistent with Ward Identities. Now, it is conve- 
nient to represent the action (17) of the system in the 
following way 



rf 3 ^(|^) 2 [(Vx) 2 +c(Vr/) 2 ] (18) 



2( 



d 3 x 



cos x • cos 



where we take 8 = to simplify things. We are looking 
for a DW solution which satisfies the following boundary 
conditions: 



(X 0,7}' 

(X -> 7T, 77' 



0) as z — \ 
■ 2tt) as z 



-00 
+00, 



(19) 



One can explicitly see from (18) that the vacuum energy 
for (x = 7r, 77' = 2-7r) at z — +00 is identically coincide 
with vacuum energy when (x, T}') fields assume their triv- 
ial vacuum values: (x = 0, r/' = 0) at z = —00. 

As we already emphasized these states (with boundary 
conditions {x = 0, 77' = 0} and {x = n, rj = 27r} corre- 
spondingly) must be interpreted as topologically different 
but physically equivalent states. Therefore, the corre- 
sponding domain wall configurations in Euclidean space 
should be interpreted as configurations describing the 
tunnelling processes rather than real DW in Minkowski 
space-time. As wc discussed in [30], a similar domain wall 
which has an analogous interpretation is known to exist 
in QCD at large temperature in weak coupling regime 
where it can be described in terms of classical equation 
of motion. These are so-called Zn domain walls which 
separate domains characterized by a different value for 
the Polyakov loop at high temperature. As is known, see 
the review papers [38, 39] and references therein, these 
Zn domain walls interpolate between topologically dif- 
ferent but physically identical states connected by large 
gauge transformations similar to our case. 

Coming back to our case, the boundary conditions (19) 
correspond to the case when the large gauge transforma- 
tion represented (in this model) by the x held is compen- 
sated by the 77' field which couples with the gluon den- 
sity operator (6) as 77' enters the action in combination 
(6 — 77') as explained above. The corresponding equations 
of motion in this simple N = 2 case take the form 
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We are looking for a solution of the system (20) which 
satisfies the boundary conditions (19). 

It is interesting to note that our system (20) with 
boundary conditions (19) formally is identical to the ax- 
ion DW configuration analysed in ref. [32] in the limit 
when the isotopical symmetry is exact, i.e. m u = rrid- In 
this case our \ field plays the role of the ir° field denoted 
as 7 in ref. [32] while \ plays the role of the axion field 
a from [32]. Precise relations are as follows: 



7T 



X 
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(21) 
(22) 



such that our equations (20) and boundary conditions 
(19) are identically coincide with the equations and the 
boundary conditions studied in ref [32]. Therefore, we 
simply formulate here the main points 2 which are rele- 
vant for our studies referring for the technical details to 
ref. [32]. 

Most important result of analysis of ref. [32] is that 
there is a unique solution of equations (20) which satis- 
fies the boundary conditions (19). The simplest way to 
convince yourself that such a solution should exist is to 
use a mechanical analogy as suggested in [32] . While ana- 
lytical formulae of the solutions is not known, its asymp- 
totical behaviour at very large distances \z\ ^> m^ 1 can 
be easily found as follows. It is clear that al large neg- 
ative z a heavy x( z ) field already assumes its vacuum 
value x = such that r)'(z) domain wall equation can be 
approximated in this region as 



VV - 2ml, sin ( ^ 



0, z< — . (23) 

m n 



Solution of this equation which vanishes at large dis- 
tances rj [z — > — 00) — > can be approximated as 



rj{z) = 8 arctan 



z <S -- 



1 



(24) 



Similarly, al large positive z a heavy x( z ) field already 
assumes its vacuum value x = tt such that rj{z) domain 
wall equation can be approximated in this region as 



V 2 t/ + 2m 2 , sin 



rf 



0, z > 



1 



(25) 



Solution of this equation which approaches its vacuum 
value 77' = 2n at large distances rj[z — > +00) — > 2ir can 
be approximated as 



rj (z) = 2tt — 8 arctan 



tan 



z > .(26) 



2 A simplest intuitive way to understand the qualitative behaviour 
of the system (20) is to use a mechanical analogy as suggested in 
[32] when variable z is replaced by time, while the fields (77', \) 
can be thought as coordinates of two particle moving in one di- 
mension with interaction determined by the potential term from 
eq.(18). 
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Formally, a similar construction when rj'(z) field inter- 
polates between different branches representing the same 
physical vacuum state was considered previously, see [34] 
for the details and earlier references on the subject 3 . 

What happens to the double layer structure (13) for 
the topological charge distribution in the presence of the 
light dynamical quark? We anticipate, without any com- 
putations, that the light dynamical field suppresses the 
topological fluctuations similar to analysis of the topo- 
logical susceptibility in this model [28]. Indeed, one can 
support this expectation by the following argument. As 
the first step we represent q(z) in form similar to cq. (13). 
The only difference in comparison with previous formula 
(13) is an emergence of the extra term due to the rf field 
on the right hand side of eq. (27), 

•m-(£)"£-W>"»(^)-<»> 

In obtaining (27) we used equations (20) and (8) to 
simplify the expression for q{z). As we already men- 
tioned, the x( z ) ne ld has a solitonic shape interpolating 
between x = and x = 71 > sec °q. (19)- Therefore, 
q(z) ~ d 2 \/dz 2 inevitably produces a double layer struc- 
ture irrespectively to the details of the solitonic shape of 
X, similar to our discussions of a pure gauge theory in 
section II C. However, the magnitude of this structure is 
strongly suppressed as a result of dynamics of the light 
quark. Indeed, the rj'(z) field assumes its central value 
rj w 7r in the region where x( z ) field varies and d 2 x/dz 2 
is sufficiently large. This factor cos \ leads to a strong 
suppression of q(z) in eq.(27) as anticipated. 

Few comments are in order. First of all, if we would 
introduce the quark in "quenched approximation" rather 
than as a dynamical degree of freedom, we would return 
to our pure glue expression (13) as the "quenched approx- 
imation" corresponds to the rf = in formula (27). Sec- 
ondly, if we would introduce a small quark's mass ^ 
into our Lagrangian it would not change qualitative pic- 
ture presented above as boundary conditions imposed on 
the system (19) can not depend on ^ 0. Indeed, 
these boundary conditions are entirely based on exact 
symmetry of the system which requires that the energy 
of the system at 8 — and 8 = 2ir is identically the same 
irrespectively to the quark's mass. Finally, the (x, r/) 
domain wall is a coherent configuration and can not eas- 
ily decay into its constituents, even though the rf is a 



3 The crucial difference with [34] is of course that the solutions 
for the system of the fields (?/,x) considered here are regular 
functions everywhere, while solution in ref. [34] had a cusp sin- 
gularity as a result of integrating out heavy fields played by the 
X field in present "deformed QCD" model. Interpretations of 
these solutions in these two cases are also very different as we 
interpret the corresponding configurations as the transitions de- 
scribing the tunnelling processes in Euclidean space-time rather 
than real static DW solution in Minkowski space-time as we men- 
tioned above. 



real physical asymptotic state of the system. This DW 
can only decay through the tunnelling process similar to 
our previous discussions, see eq. (14). Technically, it can 
be also explained using pure kinematical arguments: the 
rf constituents which are making the rf domain wall are 
off-shell, rather than on-shell, states. 

To conclude this section: the main lesson of the present 
analysis is that the deformed QCD with matter field 
supports long range correlated configurations. In differ- 
ent words, the matter and glue fields accompany each 
other in their interpolations between topologically differ- 
ent, but physically identical states. This correlation is 
enforced by very generic features of the Lagrangian (15). 
First, it is a local enforcement as the Ward Identities 
require that 8 parameter and rf field enter the effective 
Lagrangian in a specific way as cq. (15) states. Sec- 
ondly, it is a global enforcement as the 2tt periodicity in 
8 implies there existence of interpolating (X)??') configu- 
rations which inevitably present in the system. 

We emphasize once again that the long range structure 
revealed in this section might be the trace of a similar 
structure measured on the lattices, as the transition from 
high temperature weak coupling regime to low temper- 
ature strong coupling regime should be smooth without 
any phase transitions on the way [27]. However, sim- 
ilar to our comment in section II C, we expect that the 
(x, rj 1 ) domain walls determined by eqs. (20) with bound- 
ary conditions (19) and asymptotical behaviour (24, 26) 
become very crumpled and wrinkled objects with large 
number of foldings. Such fluctuations obviously increase 
the entropy of the DWs which eventually may overcome 
the intrinsic tension as holographic picture suggests, see 
[37] and references therein. In fact, it is quite likely that 
an appropriate description for this physics should be for- 
mulated in terms of holographic dual model, however we 
leave this subject for future studies. 



IV. CHIRAL MAGNETIC EFFECT (CME) AND 
OTHER TOPOLOGICAL PHENOMENA 

As we already mentioned, the ultimate goal of the 
present work is to understand the infrared physics in the 
presence of external magnetic field, which is precisely the 
environment relevant for study of the CME. We are now 
in position to couple our system (15) to external Maxwell 
U(l) field as massless quark ip carries the electric 
charge e of field. In this section we use conventional 
Minkowski metric in order to compare the obtained be- 
low formulae with known expressions written normally in 
Minkowski space-time. The rf field which appears in the 
low energy Lagrangian (15) does not couple to electro- 
magnetic field directly as it is a neutral field. However, it 
does couple via triangle anomaly, similar to the textbook 
example describing it — > 2j decay. The corresponding 
Maxwell term S~, and anomalous term Sv-rv have the 
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form 



S~ = — - / d i x 



"' 77 " 16tt 2 



(28) 



The structure of the 5'^' 77 is unambiguously fixed by 
the anomaly. It describes the interaction of the Maxwell 
field F^v with matter field. 

The interaction (28) is normally used to describe ?/ — >• 
27 decay. However, in the context of the present work we 
treat rj' (x) as external background field describing the 77' 
DW discussed in previous section III. Therefore, a num- 
ber of new and unusual coherent effects will emerge as a 
result of long range structure represented by the 77' do- 
main wall described above. Let us emphasize again that 
this long range structure being represented by the 77' com- 
ponent of the (77', x) domain wall is the Euclidean long 
range configuration describing the tunnelling processes in 
Minkowski space -time, rather than a physical configura- 
tion in real Minkowski space-time. Nevertheless, this 77' 
coherent component does interact with F^ field as eq. 
(28) dictates. 

We start by rewriting the action (28) using conven- 
tional vector notations for electric E and magnetic B 
fields, 



'77 



d^x 



1 

2 B 



Ne 2 
4^2 



i-j'E ■ B 



.(29) 



One can immediately see that in the electric field will be 
induced in the presence of external magnetic field B cx ± 
in the extended region where 77' is not vanishing, i.e. 



E 



Ne 2 
4^2 



Y-Be 



(30) 



which is precisely the starting formula in [7] if we iden- 
tify the coherent 77' component from the (t]',x) domain 
wall with induced 6 parameter introduced in [7]. Fur- 
thermore, is we assume that the DW is extended along 
(x,y) directions, we get 



( ^ ) 



(31) 



where L is size of the system, and integer number I is the 
magnetic flux of the system, J d 2 x±B* xt = $/e = 2irl/e. 
Formula (31) identically coincides with eq. (7) from [7] if 
the induced 6 parameter from that work is identified with 
extended along (x, y) configuration represented by the 
77' (z) domain wall described in section III. The induced 
electric field along z obviously implies that the current 
will flow and charge will be separated along B* xt . 



not to be confused with gluon field from eq. (4), (6) 



Anomalous coupling (29) also implies that there will 
be induced current as a result of coordinate dependence 
of the 77' field, 



Ne 2 



Stt 



(32) 



which is convenient to represent in vector notations as 
follows 

Ne 2 -> . - - Ne 2 ■ - 

Jo = — V77' • B cxt , J = — v 'B cxt . (33) 



4tt 2 



Air 2 



where we assume that the external magnetic field B cxt is 
coordinate independent. Formula (32) for the anomalous 
current has been studied previously in literature in many 
fields, including particle physics, cosmology, condensed 
matter physics. In particular, in the context of the axion 
domain wall it was extensively discussed in [32]. 

In the present context relevant for the CME formula 
(33) reduces to the well known result when the induced 
0(x,t) parameter is identified with extended 77' domain 
wall from section III. Indeed, integrating J Jod 3 x leads 
precisely to the known expression for the charge separa- 
tion effect (CSE) along z, 



Q = I d 3 xJ 



Ne 2 
4^" 



drf 
dz—— 

dz 



d 2 x ± B* 



Nle, (34) 



where we took into account the boundary conditions for 
the 7/ field (19) and replaced / d 2 x±B* xt = 3>/e = 2nl/e. 
Furthermore, the expression for J in eq. (33) can be 
presented in much more familiar way if one replaces 
rf — >• 6 = 2/i 5 as our identification suggests. With 
these replacements the expressions (33) and (34) assume 
their conventional forms which are normally used in CME 
studies. One should comment here that while our so- 
lution for the 77' domain wall considered in section III 
is time independent, in fact it actually describes tun- 
nelling effects in strong coupling regime. Therefore, it 
obviously becomes a time dependent configuration with 
a typical time scale 77' ~ Aqcd in a course of a smooth 
transition from weak coupling to strong coupling regime, 
as discussed at the end of section III. However, one can 
not use very large magnitude for 77' ~ ^5 ~ 1 GeV (as 
many people do) for numerical estimates as the effec- 
tive Lagrangian approach which leads to formulae (33), 
(34) can only be justified for small values \r/'\ <C Aqcd, 
and marginally justified for \r]'\ ~ Aqcd- For large 
7/ 1 ;§> Aqcd the effective Lagrangian approach can not 
be justified, and computations should be based on a dif- 
ferent technique when underlying QCD degrees of free- 
dom, quarks and gluons (rather than effective 77' field) 
play the dynamical role. 

One should say that there are many other interest- 
ing topological effects originated from similar anomalous 
terms as originally discussed in terms of hadronic fields 
in [7, 40], and re-derived in terms of microscopical quark 
fields in [14, 41]. In particular, if a system with chemical 
potential /i rotates with angular velocity Q, there will be 
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a current flowing along f2. The charges will be also sep- 
arated along the same direction. To get corresponding 
formulae one should replace eB — > 2/j.Q in eq. (33) as 
discussed in [7], i.e. 

Nefj, - - Ne.fi ■ - 

Jo = "TnrW ■ ^oxt , J = -TT-o-V'Qext, (35) 

which precisely coincides with eq. (A2) from [7] if one 
identifies 7/ field with induced parameter 6 from [7] as we 
already discussed above. Specific consequences of effect 
(35) relevant for heavy ion collisions were discussed quite 
recently in [42] where the effect was coined as the chiral 
vortical effect (CVE). 

The main point of this section is as follows. The long 
range structure discussed in section III might be the trace 
of a similar extended structure measured on the lattices, 
as the transition from high temperature weak coupling 
regime to low temperature strong coupling regime should 
be smooth [27]. It is important that this long range struc- 
ture describes the tunnelling effects and represented by 
(x, T)') fields in deformed QCD. These configurations are 
not real physical configurations in Minkowski space-time. 
Nevertheless, these long range configurations do interact 
with real physical ESzM field as a result of anomaly (29). 

Such an interaction transfers unphysical long range 
correlations (expressed in terms of the Euclidean config- 
urations describing the tunnelling processes) to physical 
long range correlated ESzM effects (30, 31, 33, 35). The 
corresponding coherent effects arc accumulated on large 
scales ~ L where the boundary conditions for different 
topological sectors are imposed. 

This toy model explicitly shows that the large ob- 
served intensity of the effect as studied at RHIC and 
ALICE at the LHC [1-5] might be due to a coherent 
phenomena when the large observable asymmetry is a 
result of accumulation of a small effect over large dis- 
tances ~ L 3> Aq CD . In different words, the CSE 
given by formulae (34), CME given by formula (33) and 
CVE given by formula (35) lead to a large magnitude 
for an asymmetry in spite of the fact that parameters 
\t)'\ ~ A* 5 ~ L remain paramctrically small during en- 
tire tunnelling transition. It should be contrasted with 
some other numerical estimates when the large intensity 
of the effect is achieved by choosing a relatively large 
I /ig I ~ 1 GcV S> Aqcd ~ 0.1 GeV in which case there 
is no region of validity for the effective Lagrangian ap- 
proach. 

V. CONCLUSION AND FUTURE DIRECTIONS 

The question which is addressed in the present work is 
as follows: what is the physics behind of the long range 
order which was postulated in [7], and which is appar- 
ently a required element for CME and CVE to be oper- 
ational. We attempt to answer this question using the 
"deformed QCD" as a toy model where all computations 
are under complete theoretical control as this model is 



a weakly coupled gauge theory. Still, this model has all 
the relevant crucial elements allowing us to study difficult 
and nontrivial questions which arc known to be present 
in real strongly coupled QCD. The study of these effects 
in this toy model reveals that the long range structure 
may result from the tunnelling effects and represented 
by (x,v') fields in deformed QCD. 

Apparently, such kind of transitions with long range 
structure in strong coupling regime are happening all the 
time, as it is observed in the lattice simulations. One 
should expect that the corresponding configurations at 
strong coupling should be very crumpled and wrinkled 
objects with large number of foldings in contrast with 
our smooth (x,?/) domain walls described in section III. 
Such local fluctuations are expected to occur as it pro- 
vides a large entropy for these configurations to overcome 
their intrinsic tension. If the entropy of the configura- 
tions is sufficiently large, the corresponding objects will 
have effectively vanishing tension, and an arbitrary num- 
ber of such objects can be formed and they can percolate 
through the vacuum, forming a kind of a vacuum conden- 
sate. Nevertheless, the crucial element of these DWs, the 
long range coherence, is not lost in transition from weak 
to strong coupling regime. Precisely this feature, we be- 
lieve, is the key element why the observed asymmetries 
are sufficiently strong and not washed out as a result of a 
conventional short range QCD fluctuations with a typical 
size Aqcd- 

We suspect that all other conventional mechanisms 
based on e.g. instanton/sphaleron transitions can not 
provide sufficient intensities observed at RHIC and the 
LHC as the observed asymmetries must be accumulated 
on large scales of order ~ L rather than on scales of order 
A" 1 

QCD' 

On phcnomenological side, the very basic observed fea- 
tures, such as energy and charge independence, of mea- 
sured asymmetries in heavy ion collisions are automati- 
cally and naturally satisfied within the framework based 
on long range order, see recent papers [43, 44]. The same 
framework based on the idea of a coherent accumula- 
tion of the effect also provides a natural explanation for 
a strong dependence on centrality as observed at RHIC 
and ALICE at the LHC [1-5]. 

Essentially, our study in a simplified version of QCD 
provides a precise and very specific realization of an old 
idea [45, 46] (see also [47, 48] where similar idea was for- 
mulated in different terms), that a macroscopically large 
domain with 9 in ^ ^ can be formed in heavy ion colli- 
sions. Now we can precisely identify this domain char- 
acterized by Bind ^ with interpolating long range rf 
field which traces a pure glue configuration describing 
the transition between different topological sectors. In 
different words, the domains with 6i n d 7^ should not 
be thought as real extended regions formed in Minkowski 
space-time as a result of collision. Rather it should be un- 
derstood as long range Euclidean coherent configurations 
which saturate the tunnelling transitions in path integral 
computations. Nevertheless, these long range structure 
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formulated in terms of auxiliary Euclidean configurations 
can be translated into observable long range effects (30, 
3f , 33, 35) in Minkowski space as a result of anomalous 
coupling with physical ESzM field. 

We should also mention that CME has been ex- 
tensively studied in the lattice simulations [25, 49- 
51]. Independently, very different lattice studies re- 
veal that the crucial topological configurations saturat- 
ing the path integral are represented by extended, locally 
low-dimensional sheets of topological charge embedded 
in 4d space [16-24]. Our analysis based on computa- 
tions in weakly coupled "deformed QCD" suggests that 
the long range configurations which are responsible for 
CSE/CME/CVE effects are precisely the same objects 
which we identify with long range extended objects from 
refs. [16-24]. We presented a number of arguments sug- 
gesting that this relation is in fact quite generic as it is 
based on topological features of the theory, rather than 
on a specific details of the model. Therefore, we con- 
jecture that this relation continues to hold in strongly 
coupled QCD. This conjecture can be explicitly tested 
in the lattice simulations as essentially this conjecture 
suggests that the topological charge distribution as it is 
done in [16-24] and electric charge distribution in the 
presence of the background magnetic field are strongly 
correlated and follow each other. Such a correlation also 
provides a new, and much easier way to study the origi- 
nal topological charge distribution by putting the system 
into the background magnetic field and studying the elec- 
tric charge distribution in "quenched approximation" as 
it is previously done in refs. [25, 49-51]. 

Our final comment is as follows. The transition from 
weakly coupled "deformed QCD" to strongly coupled 
regime should be smooth. Still, this transition is beyond 



the analytical control within QFT framework. What 
would be an appropriate tool to study this physics in 
strongly coupled regime? It is very possible that the de- 
scription in terms of the holographic dual model may 
provide the required tools and technique. In fact, the 
long range structure is obviously present in holographic 
model as one can see from computations of the so-called 
"Topological Casimir Effect" [37, 52] when no massless 
degrees of freedom are present in the system, but depen- 
dence of physical observables exhibit a power like sen- 
sitivity to size of the system ~ L~ p . This scaling is 
in huge contrast with ~ exp(— L) dependence which is 
normally expected in a conventional theory with a gap. 
Some analogies presented in [37] are actually suggest- 
ing that the ground state of QCD behaves very similarly 
to some condensed matter systems which are known to 
lie in topological phases. The last word whether these 
analogies can be extended to the strongly coupled four 
dimensional QCD remains, of course, the prerogative of 
the direct lattice computations. 
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